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Abstract 


The reliability of performing measurements of cosmic ray energy spectra 
with a thin ionization calorimeter has been investigated. Two Monte 
Carlo simulations were used to determine whether energy response fluctuations 
would cause measured spectra to be different from the primary spectra. 

First, Gaussian distributions were assumed for the calorimeter energy 
resolutions which (a) did not change with energy, (b) increased with 
energy, and (c) decreased with energy. The second method employed a 
detailed Monte Carlo simulation of cascades from an isotropic flux of 
protons, with incident energies representing (a) a simple power law, 

(b) a power law with a bend, and (c) an energy dependent spectral exponent. 
The results show that as long as the energy resolution does not change 
significantly with energy, the spectral indices can be reliably determined 
even for Og/E = 50%. However, if the energy resolution is strongly energy 
dependent, the measured spectra do not reproduce the true spectra. Energy 
resolutions greatly improving with energy result in measured spectra 
that are too steep, while resolutions getting much worse with energy 
cause the measured spectra to be too flat. Since thin calorimeters have 
energy resolution approximately constant with energy, they offer a viable 
means for measuring power law spectra. They are also suitable for detecting 
spectral bends or energy dependent spectral exponents, provided sufficient 
exposure is Available. 

1 . Introduction 

The energy distributions of cosmic rays carry information about 
the acceleration processes, about the sources of the particles, and about 
their propagation through and interaction with the interstellar medium. 
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As such they contain important astrophysical information about the origin 

and life history of these particles, 

X“8 

Existing data on primary cosmic rays indicate that they typically 
have energy spectra of the form % 


dN 

dE 


= K E 


( 1 ) 


Different particle species are represented by different values of K and y. 
For example, protons have a spectral index y 2,7 while iron nuclei 
may have a spectral index as low as y ^ 2,1, Examples of recent spectral 
measurements for protons and helium are shown in Fig, la, and for other 
nuclei in Figo lb. We know from previous experiments that for all species 
there is a rapid decrease in intensity with energy. Consequently, the 
major problem in designing instruments for high energy spectral measurements 
is how to achieve the maximum collecting area with a detector which has 
energy resolution sufficient for a reliable measurement of a power law 
primary spectrum. In addition, measurements should be capable of reflecting 
a spectrum which has different indices in different energy ranges; i,e. , 


dN 

dE 


^1 

E , E^ < E ^ E^ 


= E "^2^ < E ^ E^. 



Such a knee in the proton spectrum has been reported b)^ Grigprpv tp occur 

9 

just above 1000 GeV, 


It also may be that y = y(E) in some energy range, in which case 
the simple power law interpretation is invalid. We do not expect sharply 
structured features, particularly monoeBergetic lines, in the energy 
spectra of charged particles, even though they may be present at the 
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sourceo The stochastic scattering due to the propagation processes would 
spread the energies and at the same time cause the arrival directions 

to be distributed isotropically. Furthermore, we do not expect narrow- 

I 

band time structure such as the pulsations observed in radio, optical, 
and X-ray observations o The interstellar magnetic fields would spread 
the arrival times at earth of a synchronous pulse of charged particles 
emitted with small directional dispersion at the source. 

To reiterate, we expect to measure relatively smooth spectra, in 
both energy and time, for each particle species, and we expvect the energy 
spectra to be generally of the form given in Eq, 1, However we must be 
able to detect two or more. spectral indices in different energy ranges 
as indicated by Eq. 2 or even be able to detect an energy dependent spectral 
exponent. ' 

It is our objective here to show that a thin ionization calorimeter 
will suffice to measure the cosmic ray spectra in either of these cases, 
provided ample exposure time is available to cover the requisite energy 
range for distinguishing among them. Since an ionization calorimeter 
can measure primary energies over an essentially unlimited energy range, 
it is perhaps the only device which can accomplish this task in a single 
experiment. Most other energy measurement techniques are valid over 
limited energy ranges. Therefore, their use would require either several 
experiments operating in different ranges or several techniques valid 
for different ranges in a single experiment. An "ionization calorimeter" 
basically consists of an amount of matter where high energy particles 
suffer successive interactions thereby building up a shower or cascade 


origd\jai; pagp tq 

® POOR QpS 



-4- 


of particles o This way the particles transfer a significant fraction 
of their original kinetic energy to ionization of the material which 
can be measured. Thus the energy of the primary particle is determined 
by a method reminiscent of calorimetric measurements. 

Essentially, a thin calorimeter measures the energy transferred 
to electron-^photon cascades in the first interaction of the primary 
particle. The energy resolution, therefore, depends on the fluctuations 
in the amount of energy going into the production of neutral pions. The 
energ}' resolution improves as the calorimeter is made deeper because additional 
interactions occur. The additional interactions result in a larger portion 
of the incident energy appearing in the electromagnetic component so that 
at larger fraction of the incident energy is measured and thus fluctuations 
are decreased. 

Detailed Investigations using accelerator calibration, and Monte 

Carlo calculations bave been carried out on the energy response of ionization 

10-23 

calorimeters to monoenergetic particles. From that work we have 

a reasonable understanding of how the energy resolution changes with 
calorimeter depths, different absorber materials, different particle 
species, and different primary energies. A brief recapitulation of some 
pertinent results is given in Sec, 2. However, in this paper we primarily 
address the problem of how the energy resolution affects the spectral 
indices to be measured. 

We present in Seco 3 some simple analytical arguments related to 
the ptoblemo Essentially the same results were given in Refo 8. They 

p 

are repeated here for completeness and for setting the framework of' the 
subsequent Monte Carlo results using several different energy resolution 
distributions o In Seco 4p we present results of a Monte Carlo study 



using Gaussian distributions for the energy resolutionso The distributions 
were characterized by the standard deviations where is the mean 

of the distribution. We hereafter abbreviate a/iJ, as a* Calculations 
were carried out for cases (a) where q remained constant with energy, 

(b) where cr increased with energy, and (c) where ct decreased with energy. 

IiT Sec. 5. we describe results of a more direct simulation of an experiment; 

1. e., Monte Carlo-calculations were carried out for an isotropic flux 
of protons incident on a thin calorimeter whose energies were selected 
from input energy spectra. Results are presented (a) for a simple power 

law, (b) for a power law with a bend, and (c) for an energy dependent spectral 
exponent. In this latter method no specific assumptions were made about 
how the energy resolution changes with energy. The resolution for each 
energy was that predicted y- the Monte-Carlo model used to calculate 
the cascades. As pointed out in Ref, 23 the inodel is in good agreement 
with both accelerator and cosmic ray measurements of the cascade development. 
Consequently, the results given in Sec. 5. should predict the situation 
encountered in an experiment reasonably well. 

2. Energy Resolution of Calorimeters 

Investigations of the response of an ionization calorimeter to 

monoenergetic particles have shown that when a hi.gh energy hadron enters 

a large, dense absorber, almost all of the hadron’s primary energy 

is dissipated via ionization and excitation of the absorbing material. 

The area under the curve of ionization energy versus depth in the absorber 

provides a measure of E^<, For a very deep absorber the energy deposited 

can be represented by Gaussian distributions whose widths decrease as 

'*7 21 

the beam energy iricreases.'*' ~ An infinitely deep calorimeter will provide 
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energy resolution limited only by the statistical nature of the cascade 

process itself and by the frequency of sampling of the ionization energy. 

We will assume throughout this paper that the sampling frequency is 

sufficient to have a negligible effect on the resolution. The minimum 

depth which can be used for cosmic ray measurements depends on the ^rrergy 

resolution acceptable for a particular experiment. The minimum depth 

for spectral measurements would require (i) sufficient absorber depth 

so that the primary particle interacts at least once and (ii) enough 

absorber 10 radiation lengths) beyond the interaction so that ensuing 

electromagnetic cascades develop sufficiently well so that their energy 

content can be measured. The energy resolution of a thin calorimeter 

defined in this way will depend on fluctuations in the energy transferred 

(inelasticity) to neutral pions in the interaction. Our use of the L arm 

’’thin calorimeter" is not very stringent. The calorimeter used in Ref. 

8 approaches the minimum possible depth and is strictly classified as 

1-3 

thin. In previous cosmic ray exposures, the calorimeters are thin 

in comparison to the ones used at accelerators. 

The intrinsic resolution of a deep calorimeter is limited by fluctuations 

in the fraction of the primary energy going into nuclear disintegrations 

E,. of the target nuclei. As illustrated in Fig. 2 for a 3.5 \(proton 
dxs 

interaction lengths) calorimeter, this fraction decreases with increasing 
energy. For a thin calorimeter the amount of leaking out the bottom 

of the absorber plays a crucial role. The example in Fig, 2 indicates 
that the fractional energy leakage is more the higher the primary energy. 

The total unmeasured energy is the sum of the disintegration energy and 
energy leaking out the bottom. It can be seen from the top curve in Fig. 


2 that the total unmeasured energy does not drastically change with energy. 

Since the energy resolution depends among other things on the total 

unmeasured energy, the energy resolution for a thin calorimeter will 

remain approximately constant over two or more deciides in energy. This 

is illustrated in Fig. 3 where the dependence of cj oh primary energy 

is shown for two calorimeter depths. This figure represents accelerator 

2 20 

measurements with a tungsten calorimeter of about 1000 g/cm total depth. 

Notice the similarity between the shapes of the curve for the half-depth 
energy resolution from Fig. 3 and the curve for the total unmeasured 
energy in Fig, 2. The same similarity exists between the full-depth 
energy resolution from Fig. 3 and the curve for the nuclear disintegration 
energy in Fig, 2. 

The data shown in both Figs. 2 and 3 are, for events in which the 
first interaction occurred near the top of the calorimeter. In Fig, 

4 is shown a comparison of the resolutions for all incident events, and 

the subset of events for which the first interaction occurred near the 

tbp of the calorimeter. The resolutions are given as a function of the B 

calorimeter depth t for 300 GeV primary protons. It is seen that for 

both curves >he resolution improves with increasing depth, but this improvement 

is less for particles which interacted near the top. If the depth is 

sufficiently large, there is not much difference in the resolution for 

the two cases. With a deep calorimeter it matters very little where 

the first in occurs. However, for a thin calorimeter reasonable 

energy resolution can be obtained only if the particles have interacted 

and the subsequent electromagnetic cascades have developed. 

In effects a thin calorimeter measures the energy deposited by the 
first interaction of the primary particle. Therefore, its use for cosmic 

ORIGINAL PAGE IS 
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ray energy spectra measurements would result in the spectra being determined 
from a subset of the primary events. Knowledge of the interaction cross 
section and inelasticity would be required to. obtain the true flux from 
this subset. 

3 /infalytical Arguments 

The energy spectrum measured with a particular detector, in our 
case an ionization calorimeter, is a convolution of the primary spectrum 
and the energy resolution function of the apparatus. Let E represent 

f ( ' 

\r:ye primary energy of a number N(E) of events between E and E + dE, and S 
represent the represent the measured signal which is suitably normalized 
so that < S > = E for monoenergetic input. If the events of a primary 
spectrum from the interval between E and E + dE are distributed in S 
according to the di3tr''ibution f (E,S)dS, then we will call the normalized 
density function f(E,S) the "apparatus function". For a primary energy 
spectrum N(E)dE incident upon this detector, the signal spectrum actually 
registered is given by 


r 




N'(S)dS = 


(E,i3) dS dE 

E = ,0 


(3) 


We want to apply this quite general statement to the special problem 
of measuring power law spectra of th^? form given in Eq. li First we Introduce 
the new variable ' 


^ - y =■¥ 

and assume that the density function g(y) is normalized according to 


(4) 




g(y)dy = 1, with f(E.S) = gCS/E)/E 


(5) 


u 

0 
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Thus g(y) is an energy independent function which relates to the apparatus 
function in a simple way» This is not necessarily a good assumption, 
particularly not for magnetic spectrometers or Cerenkov counters where 
the spread of the apparatus function varies rapidly with energy a However, 
in case of an ionization calorimeter, where the spread is mainly due 
to the nuclear interaction properties and the inelasticity distribution, 
this is a fair assumption which is corroborated by experimental results. 

If we introduce Eqs. 4 and 5 into Eq, 3, we obtain 

OO 

1/ (S)dS = ds r N(S/y) g(y)^. (6> 

J 

y=o 


Therefore, only the primary spectrum density function N(E,), = N(S/y) still 
contains the variable S under the integral. Since the density function 
N(E) is of the form given in Eq. 1, we can rewrite Eq. 6 as 


N' (S)dS - 


dS S”'*' R 


y^'^'"^^g(y)dy. 


y=o 


(7) 


This simple results shows that the shape of the primary spectrum is reproduced 
by the detector signal, but with the normalizing constant changed by 
the factor 


M = r yW"^) g(y)dy. 

tj 

y=o 


( 8 ) 


Obviously this holds for any shape of the apparatus function f (E,S) as 
long as Eqs. 4 and 5 are satisfied and provided < S > = E. 

In general, calorimeter signals are not directly proportional to 
the primary energy. As discussed in more detail in Sec, 5, calorimeter 
signals S are more accurately given by 
< S > = aE P ■■ 


0 
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where of and P depend on the calorimeter depth but are independent of 
energy for any particular calorimeter. Following the above arguments, 
using Eq. 9 instead^, of < S > = E, we obtain for the signal spectrum 

+1) p ( V-1 ) 

N' (S)dS = ^ KS P dS y P g(y)dy. (10) 

P d 

y = 0 

This is again a power law, but with the exponent of the integral sp^ptrum 

/// 

modified by the factor 1/p, However, in practice one cap avoid this's 
situation by converting the calorimeter signal to energy using inversion 
of Eq, 9. If a and P are known from accelerator calibrations and/or 
Monte Carlo calculations, then the spectrum obtained from the energy 
measurement E , where 

c 1 ■ 

E = s (11) 

m O' 

will reproduce the primary spectrum. This case is., in fact, exactly 
analogous to the discussion leading to Eq. 7, Therefore, if the primary 
spectrum is a power law and the apparatus function does not change with 
energy, then the signal spectrum will also be a power law with the same 
exponent but with a different normalization. We should note that if y ^ 
the modified normalization integral M (Eq. 8) has a finite value for 
physically reasonable density functions g(y), whereas in case of y ^ 
the integral may or may not be divergent, depending on the exact functional 
behavior of the density function g(y) near y = 0. In practice, of course, 
we do not expect an energy spectrum with y < 1. 

4 o Gaussian Distributions for Energy Resolution 

The above analytical arguments are applicable only if the energy 
resolution distributions do not change with energy. In this section 

■' !! 
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ve present results from a Monte-Carlo study of the effect of energy 
resolutions on measurements of spectral indices. 

For simp3.icity we have chosen Gaussian distributions for the energy 
resolution. This is a reasonable approximation for calorimeters of 
moderate thickness. In addition to using distributions whose widths 
did not jfhange with energy, we considered separately cases where the 
widths increased and decreased. 

The method consisted basically of calculating-^he number of particles 

jf 

in different energy bins over two decades in energy, using five bins / 

I 

per decade. The exact number of particles in each bin was calculated 
using the chosen differential power law spectrum. Statistical f u'ctuations 
were applied to determine the number of particles in each" This 
procedure accounts for the fact that the spectrum being observSii at any 
one time is actually a statistical sample of the true primary spe:ctrum. 
Finally, the existing spectrum was convoluted with the energy resolution 
function of the detector. This was accomplished by randomly selecting ' 
the energy for each observed particle from the distribution function 
representing the energy resolution. The mean of the distribution was 

taken to be the mean energy of the particles in the energy bin containing 

i; 

the observed particle. Using the( mean energy in the bin instead of the 
actual energy of the particle reduced the computer storage and calculation 

time without significantly affecting the results. 

[■ 

Calculations were carried out for values of the spectral index y 
ij^ranging flcom 2.0 to 2.8o The measured differential spectral indices 

!i r'- ' 

/^ppear tf[-;fall in this range. For each index several widths of the energy 
resolution distributions have been considered. The distributions were 
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chosen to have standard deviations which cover the range of energy resolutions 
obtained from actual experiments. The effect of a low energy cutoff 
(e.g. geomagnetic) was tested by successively omitting low energy bins 

in carrying out least squares fits to the convoluted spectra. 

i'j 

In Fig, 5, 'are shown the results obtained fom the case where 
cr (= 0,3, 0,4, a'^^d 0,5) remains constant with energy. Each graph shows 

,;,,the value of the measured (calculated) spectral index y as a function 

..^1 u '/ 

"of the true index Yj-j-yg* solid line of linjJ.t slope indicates equal 

; O '. 

values of y and . The different sets of plotted points show 

'meas 'true r r 

the results of weighted least squares fits to the convoluted spectra. 

The open circles are for fits to all energy bins (Lq =1), while the X’s 

and the solid circles represent fits with the lower energy bins omitted 

(Lq = 2 and 3), The graphs show that for all three resolutions the value 

of Y agrees with y fot the entire range of indices calculated, 

'meas 'true ^ ’ 

provided the two lowest energy bins are omitted. For each spectral index 
value we have performed a single calculation representing a single experiment. 
The deviation of the solid points from the solid line represents the 
error one can expect in performing a single measurement of a primary 
spectrum. , ^ 

It should be emphasized that events in the individual energy bins 
of the convoluted spectrum do not represent the actual number of primary 
events in those bins. As shown in Ref. 8, the major contribl^iion to 
a particular bin stems from events of somewhat lower primary energies. 

This results from particles suffering large energy loss, and therefore 
falling into a higher bin where there are fewer particles because of 
the rapidly decreasing flux at higher energies. Since the contribution 
to a particular bin from higher energies is small in comparison to the 
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contribution from lower energies, the measured flux is higher than the 
true flux. However, there is an equilibrium in the relative number of 
particles having apparent energies higher or lower than the true energy o 
This equilibrium causes the convoluted flux to have the same relative 
energy dependence as the primary flux, i.e., the same spectral index, 
as long as boundary effects are excluded. This is illustrated in Fig. 

-6 where the fluxes are given as a function of energy for an input spectrum 
index y = 2.6 and an energy resolution <j = 0.4. The sharp low energy cutoff 
of the input spectrum destroyed the equilibrium, and consequently the 
convoluted flux is distorted near the low energy cutoff. For this reason 
it was necessary to omit the lowest energy bins in Fig. 5 in order to 
obtain the correct values for the measured spectral indices. The number 
of bins that have to be omitted depends on both the energy resolution 
and the index of the primary spectrum. 

Slowly developing cascades, with the correspondingly large amount 

of energy leaking out the bottom of a thin calorimeter, will result in 

15 23 

a much lower than average energy deposition. * This causes the energy 
resolution distributions to be not strictly Gaussians but rather more 
like Gaussians with tails extending back toward zero. In Fig. 7 are 
shown the results for the case where the energy resolution distributions 
have lead tails. Again the overall resolutions considered are 0.3, 0.4 
and 0.5, independent of energy. In order to simulate the lead tail, 
the sum of two Gaussian distributions was used. The main distributions, 
corresponding to more fully developed cascades, were used to select the 
energies of 84% of the events, while 16% of the events had energies 
selected from the wider (lead -tail) distribution. The resolution a 2 
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of the lead tail distribution was taken to be three times as large as 
the resolution of the most prominent distributiono The mean of 
the lead tail distriwas set equal to where was the mean 

of the main distribution. It can be seen from Fig. 7 that the lead tail 
does not seriously affect spectral index measurements. The main consequence 
of the lead tail is an increased disturbance in the equilibrium of the 
convoluted spectrum so that more particles are assigned energies below 
the low energy cutoff. This statement would not change if different 

« 

percentages of events had been used in the two parts of the hybrid lead- 
tail distributions. We conclude, therefore, that asymmetric energy resolution 
distributions can also be reliably used for power law spectral measurements 
provided proper caution is exercised near the cutoff energies. 

So far we have considered only distributions whose shapes do not 
change with energy. In practice, no device can measure all energies 
with the same accuracy, so the shapes will in general change with energy. 

In order to test the effect of energy dependent resolutions, we considered 
both increasing and decreasing widths of the distributions with energy. 

The results are illustrated in Fig, 8 for three monotonically changing 
energy resolution functions. The solid circles represent approximations 

2 A 

for the resolution obtained with a moderately thin calorimeter * , i.e., 

a slow increase in width as the. primary energy increases. The solid 

19 20 

triangles correspond to a deep calorimeter * , where the distribution 

widths significantly decrease with energy. The crosses represent distributions 

which become significantly wider with increasing energy. This latter 

'>} 

case is too extreme for thin calorimeters, but may be indie active of other 
types of energy measurement techniques. 
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It is clear from Fig, 8 that energy resolutions which drastically 

change with energy will result in incorrect measurements of power law 

spectral indices,, If the resolution greatly improves as energy increases 

the measured spectra will be too steep, while resolutions getting much 

worse with energy will result in spectra that are too flat. The deviation 

of the measured index v from the true index v is larger for steeper 

'meas ' true “ ^ 

spectra. This deviation is also greater the more rapid the energy dependence 
of the resolution. If the resolution is a slowly varying function of 
energy, as in the case of a thin calorimeter, the effect on the spectral 
index will be small. This is indicated by the proximity of the solid 
circles to the solid line drawn in Fig. 8. 

5. Monte-Carlo Cascade Simulations for Energy Resolution 

In this section v?e present results which do not rely on the assumption 
of Gaussian distributions for the energy resolution. Instead, we have 
carried out a complete simulation of an experiment to measure spectral 
indices. This simulation consisted of using a detailed Monte-Carlo program 
to-calculate the cascade developments for an isotropic flux of protons 
whose energies were selected from a power law distribution. The 

2 4 

calculations were specifically for a 3,5X. deep iron calorimeter ’ flown 
in two balloon exposures. The Monte-Carlo events were analyzed by the 
same method used for the flight data. The Monte-^Carlo model has been 
compared with a variety of measurements using different calorimeter con- 
figurations and absorbers. The fluctuation distributions of the energy 
loss resulting from these model calculations typically have shapes like 
Haussian distributions with lead tails similar to the ones used for Fig. 7. 
The calculations have consistently been in good agreement with measurements . 
Therefore, it is expected that this simulation technique for checking 
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the spectral index measurements should also be reliable. 

Because of the large amount of computer time required to perform 
detailed calculations for high energy cascades, it was necessary to use 
a simplified procedure. Our approach was to generate a set of isotropic 
proton events whose energies were selected from a uniform distribution 
in the logarithm of the primary energy. The calculations covered the 
range 10 - 4000 GeVo A power law spectrum was simulated by weighting 
each calculated event according to the chosen input spectrum. This technique, 
of course, does not accurately represent the statistical errors associated 
with a primary spectrum but it does permit a spectrum simulation without 
calculating an enormous number of events. Furthermore, the same calculated 
events can be used to simulate any arbitrary input spectrum. 

The convolution of the input spectrum V7ith the calorimeter resolution 

function was obtained in the following ways. First, the number of 

calculated events in various energy bins was determined, using 5 bins 

per decade for the input Monte-Carlo energies E^, Next the energy 

measured by the calorimeter for each event was determined from the algorithm 

(discussed below) used to convert the measured signal to primary energy. 

The weighting factor used for each event having signal '^^as K E^ > 

where N was the number of events actually calculated for the energy 
3 

bin in which the true energy E^ was located* The convoluted flux was 
then given directly by binning this weighted flux in terns of E^, also 
using 5 bins per decade,, 

A schematic diagram of the GSPC calorimeter used in this simulation 
is shown in Fig. 9. It consisted of 7 modules, each 0.5 proton interaction 
lengths thick, with a total depth of 3,5 interaction lengths. The modular 
stiaicture easily affords testing calorimeters of different depths o By 
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neglecting signals in the downstream modules we have considered depths 
from 2.0 to 3.5 interaction lengths. 

The response of the calorimeter to particles of primary energy 
is well represented by the relation 

The Calorimeter signal is the summed number of cascade particles 

in k modules, where the summation begins with the module in which the 
first inteijaction occurred o The parameters a and g used for this work 
were ay = 6.4, =0.97, a^ =5.7 and 3^ = 0.9. These values were obtained 

from a least squares fit to a plot of log 3sl vs log for all the 2087 
cascades calculated in detail as the first step in this simulation procedure. 
They are valid for an isotropic flux of protons whose trajectories are 
restricted to pass through both the top and bottom of the calorlmetero 

Therefore, for each calorimeter depth (k modules) the estimated energy 

"k. ■ ■ 

E^ assigned each particle was determined from the algorithm 
m , 

Our results for energy spectra measurements, using the cascade simulations 
for primary protons, are given in Figs. 10, 11, and 12. Figure 10 is 


for a simple power .law, while Figs. 11 and 12 are for a power law with 




a bend and an energy dependent spectral exponent, respectively. In each 
figure the logarithm of the differential flux is shown as a function 
of primary energy E^. The circles represent the Monte-Garlo input flux, 
the crosses represent the flux obtained using 4 modules for the energy 
measurement, and the triangles represent th? flux using energy measurements 
from all 7 mdduleso Notice that the flux from 4 modules is higher than 
the input flux for energies considerably greater than the cutoff energy 
(bin 6). The reason for this effect was explained above in relation to 
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Fig. 6. The flux obtained from 7 modules' is less than the input flux 
.because of the smaller total number of events that pass the selection 
criterion: Not all events had cascades developing in all 7 modules » 

However, if suitably nor^J-’llzed, this flux would also be greater than 
the input fluxo 

It seems to be clear from Figs. 10, 11 and 12 that the calorimeter 
signals are capable of reproducing the input spectrum, whether the spectrum 
is a simple power law, a power law with a bend, or an energy dependent 
function. This is true for energy measurements using either 4 modules 
or 7 modules of the calorimeter. This result substantiates our earlier 
■'conclusion frQm Fig, 6 that a moderate energy dependence of the energy 
resolution function will not greatly distort spectral measurements. 

It should be pointed out here that the proton energy resolution for 4 
modules of the GSFC calorimeter ranges from 32% to 38% as the primary 
energy goes from 17.6 GeV to 1000 GeV. For the same energy range the 
resolution for 7 modules changes from 20% ip 23%. 

In calculating the results given in Fig. 11, we arbitrarily set 
the abrupt change in the spectral index from 2,5 to 3.0 to occur at 
500 GeV. This energy value was chosen in order to have data for several 
energy bins on either side of the kirik. Qualitatively, the results would 
not differ regardless of where the index change occurred. However, as 
a practical matter, it would be necessary in an experiment to have energy 
measurements over about one decade in energy above such a kink in order 
to detect it. 

The necessity of a large energy range is even more acute in case 
of an energy dependent spectral exponent. In fact, straight lines 
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exemplifying simple power laws, provide very good fits to any 5 points 
(one decade) of the flux presented in Fig. 12, as long as bins near the 
low energy cutoff are omitted. Experimentally, for example, one might 
interpret that energy dependent spectral exponent as a spectrum consisting 
of different power law indices on either side of a kink at about 500 
GeV. Consequently, differentiation between a spectrum with a bend and 
an energy dependent exponent requires at least three and preferably four 
decades of energy coverage. 

6. Conclusion 

Evidence has been presented to show that a thin ionization calorimeter 
can be used to determine the energy spectra of high energy cosmic rays. 

Since many primary particles may penetrate a thin absorber without interacting, 
a thin calorimeter is useful only for measuring energies of particles 
that interact hear its top. Consequently, one must know the relevant 
interaction cross sections in order to determine absolute particle fluxes. 

Since many cross sections measurements have been, or can be, made using 
the high energy accelerators now available, this does not in principle 
limit the use of thin calorimeters for energy spectra measurements. 

Regarding the effect of poorer energy resolution on measured energy 
spectra, it has been shown that the index of a single power law spectrum 
can be obtained with energy resolutions of 50% or worse. In order to 
detect structure such as kinks and bumps in the spectra it is desirable 
to have both good energy resolution and good statistics. A bump in a 
spectrum which has a width smaller than the energy resolution will appear 
smeared out in the measurements. An abrupt kink will also tend to smear 
out more as the energy resolution decreases. If a spectrum really consists 
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of two separate power laws, it is necessary that good statistics be obtained 
for events with energies substantially greater than the energy where 
the kink occurs. The coverage of a large energy range in a single experiment 
is necessary to determine whether the primary spectral exponents are 
actually energy dependent. 

In experiments involving the low fluxes of very high energy cosmic 
rays one must invariably make a compromise between the accuracy of energy 
measurments and apparatus collecting power o All calorimeters used in 

cosmic ray experiments near the top of the atmosphere have been in fact, 

1 2 4 8 8 

relatively thino * * * The one used by Schmidt et al. approaches the 

minimum allowable depth as defined in Sec. 2. There are alternate energy 

measuring techniques known that allow instruments with large geometry 

factors to be built such as gas Cerenkov detectors and transition radiation 

detectors. These techniques suffer from the fact that the dynamic 

range is relatively narrow for each individual Instrument 

and a spectrum would havp to be pieced together using several different 

Ij 

techniques. In piecin^^' data together from several instruments one can 
never rule out that syS’tematic errors enter that, while giving satis- 
factory values for the fluxes at different energies, might partially 
obscure the Information that is contained in the detailed knowledge of 
the spectral shape,, A ' bin ionization spectrometer seemsS to be the only 
currently practical means tr determining the shape of cosmic ray particle 
energy spectra over a ve *ge dynamic range in one single experiment. 

Deeper calorimeters prov 'e accurate energy measurements, but being 

subject to always present •'nd cost ] 'tations, the deeper calorimeters 

have smaller collecting arec lid angles. Our studies 
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have shown that measurements of the spectra of high energy particles 
would be more affected by the low numbers of eVents collected, thereby 
limiting the energy range covered, rather than by the accuracy of the 
energy measurements of those events. For this reason, a thin ionization 
calorimeter should provide a viable means of answering timely questions 
on the energy spectra and composition of high energy cosmic rays. 
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Figure 
Figc i 

Fig. 2 


Fig. 3 


Fig. 4 


Fig. 5 


Captions 

2 

(a) -.Energy spectra of protons and helium nucleio- 

3 

(b) Energy spectra for heavier nuclei. 

Typical results for the unmeasured energy in an ionization 

calorimeter as a function of the primary proton energy E . 

These results are from a Monte Carlo simulation of the cascade 

development for 3.51 deep iron calorimeter. 

(See Sec. 5/o Nuclear disintegration energy is 

represented by E^^^ while the energy leaking out the 

bottom is given by E ; the sum E,. + E constitutes^ 

^ out das out '"■'sv 

the total unmeasurd energy. ''' 

Energy resolutions for two calorimeter depths as a 

20 

function of primary energy E^ . The energy measured by 

the calorimeter is given by E . These results are from accelerator 

2 

measurements with a tungsten calorimeter of 1000 g/cm total depth. 
The abbreviation represents the "high resolution module" 
located at the top of the calorimeter. 

Energy resolutions as a function of tungsten calorimeter depth 

20 

for 300 GeV protons. For the upper curve the location of , .the 

first interaction could be anywhere. The lower curve includes only 

events in xr^hich the first interaction occurred in, the high 

resolution module (HRM) at the top of the calorimeter. 

Calculated results for the measured spectral indices v as a 

' meas 

function of the true indices v.. ° Gaussian distributions were 

' true 

used for the energy resolutions, (a) Relative standard deviation 
cr = 0.3, (b) cr = 0.4, arid (c)^ct = 0.5. Weighted least squares 


Q 
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fit for all energy bins is indicated by L =1, while L = 2 is 

<;• 

Q 

for fit omitting lowest energy bin, and = 3 is for fit 
„ omitting the two lowest energy bins. 

Figo 6 Typical comparison of the input power law flux with the flux 

obtained after folding in the energy resolution distributions. The 

calculated flux has 8645 events below the -low energy cutoff o 

(I .f 

Fig, 7 Calculated results for the ^measured spectral indices y a 

'Tn0ds 

function of the true indices . Two Gaussian distributions 

were used to simulate symmetrical energy resolution distributions 
having lead tails. The overall relative standard deviation a was 

taken to be energy independent. The main distribution (a^) was 

O' '■ 

used for 84% of the events, and the lead tail distribution 

(CT 2 ) of events. 

Fig. 8 Calculated^; results for the measured ‘spectral indices' y as a 

function of the true indices . These results are for 

' true 

Gaussian energy resolution distributions whi,ch have the indicated 

energy dependence. The solid line represents equal values of 

V and V o The dashed lines are visual fits to the data, 
•meas 'true 

Fig, 9 Schematic diagram of the iron calorimeter used by the Goddard 
' Space Flight Center. The results in Figs. 10, 11, and 12 were 
calculated for this calorimeter design. 

Figo 10 Calculated results for the differentiaT fluxes measured with an 

iron calorimeter of two different depths. The input spectrum, was 
a simple power-law with spectral index y =2.7. The measured 
fluxes were determined from the cascade':^energy deposited in 4 
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modules and 7 modules, respectively, of the iron calorimeter 
illustrated in Fig. 9. Parallel lines were drawn through the 
calculated points to show that the measured spectral indices, 
represented by the slopes of the lines, reproduced the input 
spectral index. 

Fig. 11 This plot is the same as Fig. 10, except the input spectrum index 

* ; 

was changed from 2.5 to 3.0 at 500 GeV. 

Fig. 12 This plot is the same as Fig. 10, except the input spectrum had an 
energy dependent exponent. Note that, although the data points 
representing measurements follow the curvature of the input data, 
straight lines also provide good fits over a single decade (five 
data points) in energy. Over two decades in energy, two straight 
lines of different slopes, analogous to Fig. 11, would be a good 
fito 
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MONTE CARLO PRIMARY ENERGY 
ENERGY FROM SUM OF 4 MODULES 
ENERGY FROM SUM OF 7 MODULES 
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-O- monte carlo primary energy 

ENERGY FROM SUM OF 4 MODULES 
\ ENERGY FROM SUM OF 7 MODULES 





o MONTE CARLO PRIMARY ENERGY 
X ENERGY FROM SUM OF 4 MODULES 
A ENERGY FROM SUM OF 7 MODULES 


y = 2.25 + 0.25 LOG E 






